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Abstract 


Assumption  of  the  uniform  gravity  field  was  justified  by 
computing  errors  due  to  this  assumption  and  establishing  the 
limits  within  which  the  simplified  computations  should  be 
performed.  Vertical  take-off  and  ascent  was  analyzed  for 
different  values  of  a  constant  propellant  flow.  The  analysis 
of  the  vertical  descent  and  landing  was  extended  to  the  cases 
of  a  constant  propellant  flow  and  thrust,  to  the  modulated  and 
intermittent  thrust.  Typical  numerical  examples  were  computed 
without  embarking  on  a  generalized  optimization  problem,  but 
indicating  the  ways  of  minimizing  the  propellant  consumption. 


iii 


CHAPTER  I 


Physical  Conditions 


1.1  The  Moon  is  not  an  exactly  spherical  body.  It  has  an 
equatorial  bulge  pointing  towards  the  Earth,  hence  its  shape  can 
be  approximated  as  an  tri-axial  ellipsoid. 

In  the  present  report.,  however,  the  Moon  is  considered  as 
a  spherical  body  of  a  radius  R.  There  is  no  sea-level  on  the 
Moon  and  lunar  ground  level  gravity  acceleration  corresponds  to 
the  level  of  large  maria,  like  Mare  Imbrium.  The  slow  axial 
rotation  of  the  Moon  allows  to  neglect  the  centrifugal  effect  on 
the  gravity.  The  following  data  have  been  accepted  for 
computations : 

Earth's  true  equatorial  gravity  ge'  =  32.2284  ft/sec2 


Earth's  standard  gravity  measured  „ 

at  \  =  46°  g  =  32.174  ft/sec^ 

e 

The  mass  of  the  Moon  Mm  =  0.0123  Me 

Lunar  ground  gravity  acceleration  g0  =  0.164  ge' 

o 

Lunar  ground  gravity  acceleration  g0  =  5.3  ft/sec 

Lunar  mean  radius  R  =  1080  stat.  miles  =  1740  km. 
Sidereal  period  of  axial  rotation  T  =  27.32166  days 


Angular  velocity  of  rotation  co  = 


2t r 
T 


.  _ n  rad  .  .  _  . 6 

0.23  -  2.66  x  10 

day 


rad 

sec 


Tow  velocity  at  the  equator  o>R  =  2.88  x  10 


-3  mi les 


15.2 


feet 


sec . 


sec . 
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2  ^ 

Gravity  constant  g0R  =  kM 


1.17  x  103 


(mi  les)3 
- 9 — 

sec* 


Square  of  ground  circular  velocity  V 
Lunar  ground  circular  velocity  Vc  = 
Lunar  ground  escape  velocity  Ve$  = 


c2  =  g0R  =  1.084  2^ 

1.041  mi/sec  =  5500  ft/sec 
1.472  mi/sec  =  7770  ft/sec 
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CHAPTER  II 


Errors  Due  To  Uniform  Field  Assumption 


2.1  Assumption  of  a  uniform  gravity  field  for  the  final  phase  of 
the  lunar  landing  or  for  the  initial  phase  of  the  take-off,  can 
be  justified  if  we  establish  the  volume  within  which  the  error 
due  to  this  assumption  is  negligible.  If  wc  extend  the  simplifying 
assumptions  beyond  this  volume  a  method  should  be  provided  to  find 
easily  the  corrections  for  computed  values  of  range,  velocity  etc. 

Let  us  examine  first  the  variation  of  the  lunar  gravity 
acceleration  in  the  vertical  direction.  As  there  is  no  sea-level 
value  for  this  acceleration,  what  we  call  a  ground  value  gQ 

may  not  have  the  same  value  on  the  lunar  surface.  In  the  present 
paper  these  local  divergencies  are  neglected  and  the  ground  level 
of  g  is  assumed  to  be  g0  =  5.3  ft/sec2  .  Assuming  sphericity  of 
Ahe  Moon  and  concentric  density  distribution  we  express  the  gravity 
variation  with  the  altitude  y  over  a  hypothetical  ground  level. 


9  =  «0 


=  g, 


2. 3. 4. 5 


1 

R 


Assume  as  a  mean  lunar  radius  R 


1.2. 3. 4 

1080  stat.  miles  and  find  the 


maximum  relative  error 


42 

g0 


in  replacing  the  variable  g  by  a  constant 


value  g0  up  to  the  altitude  y  =  10.8  st.  miles  =  57,000  ft. 
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Ag  _  flo  ~  g  _  2  l 


3  y 

1  -  -  -|+2 
21  R 


5 

2 


if  we  take  ^  =  - —  , 

R  100 


Ag  _  2  _ 

g0  100  l  2  x  102  104  2  x  106 


3 

1  -  — 


+  .  • 


=  I__  I  1-Q.015  +  0,002  -  0.0000025  +... 
100 


-2-  x  0.9852  =  0.0197 
100 


hence  the  error  is  slightly  less  than  2  percent  of  the  initial 
value.  By  accepting  this  error  we  should  investigate  the  time  and 
altitude  errors  due  to  this  approximation. 


2.2  Now  let  us  consider  the  errors  in  horizontal  motion  (range 
error )  . 

In  the  uniform  gravity  field  we  assume  also  that  a  planet  has 
an  infinite  radius,  and  hence  its  surface  is  an  infinite  plane 
passing  through  the  point  selected  as  an  origin.  The  sphericity  of 
the  planet  produces  the  inclination  of  local  vertical  with  respect 
to  OY  at  the  origin.  The  range  is  measured  along  the  plane  tangent 
to  the  surface  at  0,  and  the  curvature  of  the  surface  may  increase 
it  as  shown  on  Figure  1.  The  point  P  where  a  landing  may 
theoretically  occur  is  actually  at  an  altitude  Ay  above  the  ground 
level.  For  any  arbitrary  horizontal  distance  x  let  us  compute 
the  angular  (Ay)  and  altitude  (Ay)  errors. 
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tan  Ay  =  - 
R 


or 


2  2  2 
<R  +  Ay)  =  R  +  x 

2  =  xl 

R  R2  R2 


neglecting  the  higher  powers  of 

Ay  _ 

R  2R2 


R 


' 

we  finally  have 

i 


TABLE  I 


x 

2 

4 

6 

8 

10 

12 

14 

Stat.Ml . 

5  *  103 

1.85 

3.70 

5.56 

7.41 

9.26 

11.11 

12.97 

ay 

0°06 ' 

0°  13* 

0°19’ 

0°25 ' 

0°32 * 

0°38 ' 

0°45  * 

Ay  .6 

-  x  10 

R 

1.71 

6.85 

15.45 

27.45 

42.87 

61.72 

84.11 

Ay  x  103 

1 .848 

7.4 

16.69 

29.62 

46.28 

66.62 

90.90 

St .Mi  . 

Ay 

9.8 

39.1 

88.1 

156.5 

24  4 

352 

480 

Ft. 

From  the  Figure  2  an  approximate  extension  of  the  range  due  to 
sphericity  can  be  estimated  a»  Ax  =  Ay  cot  «C,  where 
«c  =  180°  -  *<F, 


«<F  is  the  fall  angle. 
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Figure  2 


I 
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2.3  The  errors  in  the  estimation  of  time  and  velocity  are  computed 
only  for  the  case  of  a  free  vertical  fall. 


or 


In  a  uniform  gravity  field 

y  =  ”  0o  J 


d  y  ,  dy 

— r,  =  y  —  =  -  g0  . 

d  t  *■  d  y 


The  integration  yields 


ii-io! . 


-  o„  <J  -  r„>  . 


If  the  initial  speed  y0  =  0,  and  if  y  =  0  that  is  the  fall 
ends  at  ground  level  we  obtain  the  impact  velocity 

*i2  =  2  goyo  * 

Consider  now  a  free  fall  in  a  variable  gravity  field,  obeying 
inverse  square  law. 

.  dy  R2 

y  —  “  -  90  - o  ’ 

dy  (R  +  y)2 


hence 


y2  -  y02  2 

- ^ -  =  g0R 


1 


R  +  y  R  +  y0 
If  y  =  0  and  y0  =  0,  we  obtain 


?2  =  2  -**- 


=  9oR 


2  y0  -  y 


(R  +  y)(R  +  yfl) 


R  +  y0 

The  difference  between  the  vertical  impact  speed  yy  of  variable 
field  and  that  of  the  uniform  one  will  be  found  from  the  equation 
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R  2  a  v  ^ 

y2  -  yv2  =  (y  -  yvHy  +  yv  )  =  2  g0y0  l - -  - 

L  R  ♦  y0  J  R  +  y0 

Hence  the  error  due  to  simplified  assumption  is  of  the  order 


By  using  serial  expansion  the  error  function  can  be  given 
the  following  form 
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For  the  Moon  -yj2  gQR  =  escape  velocity  from  the  ground  level 
=  7770  ft/sec. 

The  following  table  of  impact  velocity  errors  has  been  prepared 
assuming  a  free  fall  from  the  altitude  y0  with  the  initial  zero 

velocity.  The  correct  values  of  impact  vertical  velocity  have 
also  been  computed  from  the  formula 


and  the  relative  error  estimated  ... 


TABLE  II 


8.1 

10.8 

16.2 

21.6 

42800 

57100 

85500 

114000 

.0075 

.0100 

.015 

.02 

1.00374 

1 .00499 

1.00747 

1.00995 

3.77 

5.03 

7.56 

10.1 

670 

774 

945 

1090 

x  10 


Ft/sec 
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So  the  velocity  error  is  usually  less  than  one  percent  for 

y0  <  100,000  ft.  If  slide  rule  is  used  the  difference  Ay  should 

be  estimated  from  the  formula  and  not  from  the  difference  y  -  yv  . 

In  the  latter  case  the  subtraction  of  two  almost  equal  numbers 
greatly  magnifies  the  error. 


2.4  The  time  of  free  fall  from  the  altitude  y0  in  the  uniform 
gravity  field  is  estimated  from  the  formula 


in  seconds  if  the  initial  zero  vertical  velocity  is  assumed,  and 
if  y0  is  given  in  feet. 

In  the  variable  gravity  field  the  time  tv  is  estimated  from 

the  Kepler’s  equation  applied  to  a  rectilinear  ellipse,  that  is 
with  an  eccentricity  e  approaching  unity:  e— >1.0.  An  assumption 
is  made  that  in  the  vicinity  of  ground  level  the  velocity  is  less 
than  its  escape  value,  then  t  and  E  are  measured  from  the  apofocus 
of  an  elliptic  orbit 


n  I tv  -  t0 


E  +  e  sin  E, 


where 


n2  = 


and 


the  semi-major  axis  a  can  be  defined  from  the  initial  conditions 

g  R2 

0 


2  g0R 


R  +  y, 


-  v. 
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and  as  V0  =  0,  we  obtain  a 

at  the  Moon's  center. 

8  90r2 


R  +  y, 


that  is  the  perifocus  lies 


n2  = 


n  = 


R  + 


yo 


vr 


yo 

1  4  — 


3/2 


Substitution  of  numerical  values  yields 


Yj> 

R 


-3/2 


•3/2 


1 

sec 


For  e  =  1  the  ellipse  is  reduced  to  a  straight  line,  the 
eccentric  anomaly  E  is  computed  from  the  Figure  4. 
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At  the  ground  level 
R  +  y 


ON*  =  a  -  y  = 
3  o 


R  -  y. 


2 

R  -  y0 


-  y«  = 


ON 

cos  E  =  —  = 

a  R  +  y0 


sin  E 


1  -  cos  E  = 


E  =  arc  sin 


-V 

_  2  V  n/« 

~  - . -  t 

i  + 

2  V  y0/R 


i  - 


R  -  y< 

R  +  yf 


*  yyj 

»  +  yn 


1  *  y./"- 

Substitute  into  the  Kepler’s  equation,  taking  tQ  =  0,  and  e  =  1, 


tv  =  -  I  E  +  e  sin  E  I  =  3.667  x  10“  x 


1  +  12- 
R 


3/2 


E  +  sin  E 


The  error  due  to  the  uniform  field  assumption  is  given  by 
the  difference. 

3/2 


At  *  iv  -  t  s  -  0.614  3.667  x  102  x 


1  + 


yo 


E  +  sin  E 


The  results  are  given  in  the  Table  III. 
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At/t  0  3.85  7.8  8.1  19  x  10“3  (14) 
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2.5  Having  established  the  limits  of  errors  due  to  the 

assumption  of  a  uniform  gravity  field  we  may  consider  a  cylindrical 
box  whose  base  is  resting  on  the  lunar  surface.  The  initial  point 
of  any  trajectory  lies  on  the  vertical  axis  of  symmetry  of  the 
box;  whose  radius  is  14  stat.  miles  =  74000  ft.,  its  height  =  22 
stat.  miles  =  114000  ft.  If  the  initial  and  end  point  of  a 
trajectory  lie  within  this  box  the  errors  due  to  simplifying 
assumptions  are  limited  by  the  following  values: 

Ay  <  0°45 ' 

Ay  <  480  feet  (altitude  of  the  impactfor  horizontal  range) 


Ay 

y 

At 

t 


< 


< 


1  x  10"2 

2  x  10“2 


If  a  further  sacrifice  of  accuracy  is  permitted  by  extending  the 
dimensions  of  the  box,  the  following  formulas  enable  us  to 
estimate  the  order  of  magnitude  of  the  error. 

Ay  _  x^ 

R  ~  2  r2  ’ 


Ax  =  (Ay)  cot  (180°  -=  <<p) 


=  3885 


1 12  ' 
R 


3/2 


ft/sec . 


3/2 


E  +  sin  12  -  0.614 


TR  63-2,  Pg.  17 


CHAPTER  III 


Powered  Ascent 

3.1  The  total  thrust  of  a  rocket  engine  is  given  by  the 

expression 

F  =  "  u  +  Ae  |  pe  -  p0j  (3-1) 

where  m  is  the  mass  flow  of  the  propellant  in  slugs  per  sec 
u  the  actual  exhaust  velocity  in  ft/sec  measured 
relative  to  the  rocket  at  the  exit 

n 

Ae  =  the  exit  area  in  ft6 

p  =  the  exit  pressure  in  lb/ft6. 

© 

p0  =  surrounding  pressure 

In  the  vacuum  p0  =  0  and  the  thrust  assumes  a  value 

F  =  m  u  +  A0  pe  (3-2) 

which  can  be  replaced  by 

F  =  uef f . 

where  ueff  is  the  effective  exhaust  velocity  found  from  the  equation 

^  A  e  p  g 

u  -,  =  u  +  — -  -  (3-3) 

ef  f  • 

m 

In  further  development  the  index  eff.  will  be  dropped  and  u  will 
always  mean  an  effective  velocity.  Hence  F  =  m  u. 
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The  expression  for  the  specific  thrust  I  of  a  rocket  is 
obtained  from  the  equation: 

F=mu=Iw=mgeI  ;  (3-4) 


where : 


F  —  the  total  thrust  of  the  rocket  engine  in  pounds. 

u  —  the  effective  exhaust  velocity  at  the  nozzle's 
exit  in  ft/sec.  The  term  m  u  includes  the  exit 
pressure  term. 

m  —  the  mass  flow  of  the  propellant  = 
fuel  +  oxydizer  in  slugs/sec. 

I  —  specific  thrust  in  — — — (Thrust) ___  _  $ec< 

lb/sec  (propellant) 

w  --  the  propellant  consumption  in  lb/sec. 


hence 


»  =  I  !  =  I  =  I  g 

mm  e 


(3-5) 


ge  —  is  a  constant  coefficient  giving  the  ratio 

between  the  specific  thrust  and  exhaust  velocity 


(I  =  — ),  it  is  not  affected  by  varying  gravity 
9e 


field. 

gg  is  not  the  true  gravity  acceleration  but  is  the  standard  NACA 
gravity  acceleration  corresponding  to  that  measured  on  the 
surface  of  the  Earth  at  approximately  X  =  46°N,  hence  it 
corresponds  to  the  conditions  at  which  rocket  fuel  weight  is 
measured  on  the  Earth.  Its  value  =  ge  =  32.174  ft/sec^. 
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3.2  Consider  now  the  vertical  ascent  of  a  rocket  in  airless 
conditions.  The  forces  acting  at  every  instant  on  the  remaining 
mass  M  of  the  rocket  are  those  of  the  thrust  and  weight,  hence 


M 


d V 
dt 


F 


Mg 


In  this  problem  a  constant  g  =  g0  is  assumed,  provided  the 

altitude  is  below  an  imposed  value,  the  exhaust  effective  speed 
u  is  also  assumed  to  be  constant,  hence: 


M 


dV 

dt 


m 


u 


^00  9 


or 


dM 

u  —  -  Me  ,  as  m 
d  t  o 


dM 

dt  ’ 


hence 

...  dM 

dV  =  -  u  —  -  g  dt  ; 

M  o 

The  integration  yields: 

Mo  Mq 

V  -  Vo  =  u  In  r  -  go  t  =  I  ge  In  r  -  go  t 

If  VQ  =  0,  the  equation  assumes  the  form: 


(3-6) 


.  M 

V  —  u  In  — —  - 
M 


90  ** 


(3-7) 


where  In  is  the  natural  loge  =  2.3026  log10  ,  g0  has  been 
defined  as  the  local  ground  gravity,  M0  is  the  initial  mass  of 

the  rocket  at  time  t  =  0,  M  is  the  terminal  mass  of  the  rocket 
at  time  t  =  t.  The  ratio  of  masses  can  be  replaced  by  the  ratio 
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of  weights  whatever  the  gravity  field  is 


(3-8) 


Insofar  no  assumptions  have  been  made  as  to  the  rate  of  the 
propellant  consumption,  and  so  the  equation  (3-6)  and  (3-7)  can  be 
applied  to  the  cases  of  constant  or  variable  fuel  flow.  Only 
a  constant  u  and  consequently  a  constant  I  has  been  postulated. 

If  the  rate  of  the  fuel  flow  is  constant  and  equal  to 
a  fixed  fraction  of  the  initial  mass,  we  write 


m  =  K  MQ,  and  the  mass  M  after  t  seconds  since 


the  ignition  becomes 

M  =  (1  -  Kt ) ,  conversely  W  =  WQ  (1  -  Kt) 


3.3  Assume  now  that  the  rocket  ejects  a  constant  mass 

m  =  K  M0  per  second,  where  M0  is  the  initial  mass  of  the  rocket. 
Let  M i  be  the  mass  of  the  empty  rocket,  hence  M0  =  Mj  +  fuel. 

Show  which  conditions  must  be  fulfilled: 

a)  that  the  rocket  can  rise  at  once  after  the  ignition. 

b)  that  it  can  rise  at  all. 

If  it  can  rise  vertically  at  once  find  its  maximum  velocity 
and  maximum  height,  assuming  a  constant  gravity  acceleration  g0. 
A  rocket  can  rise  at  once  if  its  thrust  is  greater  than  the 
initial  weight 

F  >  W„ 
or 


K  M0  u  >  Mq  g0 
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9  o  9  o  v 

K  «  >  50  or  I  ge  >  -  ,  I  de  *  -  >  0 

It  can  rise  after  some  time,  if  its  thrust  is  at  least  greater 
than  its  empty  weight 


Assuming  K  u  > 
of  ignition, 
at  any  moment 


K  M0  u  >  Mj  g0 

Mj 

K  «  >  u;  «.  • 

g0,  we  can  compute  the  ascension  from  the  moment 
From  the  previous  sections,  the  vertical  speed  is 


»  Mo 

v  =  vo  "  9o  1  +  u  jj“  : 

If  V0  =  0  and  M  =  M0  -  K  M0  t  =  M0  (1  -  Kt )  , 

then 

V  =  -  g0  t  -  u  -tn  (1  -  Kt)  =  -  gfl  t  -  I  Qe  In  (1-Kt) 

(3-9) 

The  acceleration  at  any  time  t  of  powered  flight  is 


dV 

dt 


dV 

dt 


K  u 

1  -  Kt 


and  the 
second 
deri vati ve 
’  of 
V 


K2  u 


(1  -  Kt)' 


The  mathematical  extremum  of  velocity  is  attained  when 


0,  the  corresponding  value  of  t  is 


1  =  Ktm  = 


Ku 

So 


1  u 


this  value  corresponds,  however, 
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to  a  minimum  value  of  V,  as  the  second  derivative  of  V  shows:  this 

dV 

is  apparently  contradictory  to  the  fact  that  —  is  positive  and 

d  t 

increasing  with  t. 

The  paradox  is  explained  by  taking  notice  that  tm  is 
1  u  g  -  K  u 

negative,  because  -  -  — *  =  •  <  0.  So  this  minimum  V  lies 

K  g0  K  g0 


in  the  past  time  before  the  ignition  (and  motion)  started. 

During  the  ascension,  insofar  K  u  >  g0  the  velocity  and 

acceleration  (in  the  vacuum)  are  both  increasing  until  the 
exhaustion  of  fuel,  hence  the  maximum  value  of  the  velocity  will 
be  attained  at  a  certain  time  tj,  when  M0  (1  -  Ktj)  =  Mj,  therefore 

M1 

1  -  Kti  =  iji  . 


and 


1  Mi  1  Wi 

t !  =  -  ( 1 - -  )  =  -  (1 - i  ) 

K  M„  K  W„ 


Hence 


max 


5»  M  "l,  *  >  "» 

-  r  °  ■  u  tn  in 


9  o 


(1  -  — )  +  I  ge  tn  —  .  /o  10) 

K  Wo  e  Wl  (3-10) 

The  value  of  the  acceleration  at  this  moment  also  attains  its 
maximum  value: 

<£>  =  -  5  . K  M«  " 


dt 


M 


K  W, 


I  ge 


=  -  g0  + 


(3-11) 
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Hence 


(dV/dOj 


-  0.164  +  K 


I  ( 


Since  this  moment  the  positive  term  of  the  acceleration 
vanishes  and  its  value  becomes 


dV 

dt  =  "  V 


3.4  The  computation  of  a  trajectory  consists  of  two  parts: 

powered  flight  and  coasting.  Assume  the  initial  conditions 
To  =  °»  vo  =  °* 

In  powered  flights 

dy  .1 

V  =  — ■  =  -  g0t  +  u  In  - - —  ;  assuming  that  the 

at  1  -  Kt 


take  off  was  at  the  ground  level,  i.e.,  y  =  0,  we  obtain  the 
altitude  y  at  any  time  t  after  the  ignition. 


y  = 


2 


1 

( - )  dt, 

1  -  Kt 


ct  . 

f  -tn(l-  Kt)  dt 
Jo 

=  t  In  (1  -  Kt)  - 

=  t  -tn  (1  -  Kt)  - 

=  (t  -  i)  In  (1  - 
K 
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hence 


g0t2 

y  ~  2 


u  -  (t  -)  In  (1  -  Kt)  +  t  (3-12) 

K 


Let  the  altitude  at  burn-out  be  y0;  the  value  of  y  can  be  found 


by  setting 


=  t  =  i  (1  -  %  , 

1  K  M.  ’ 


M1 

1  -  Kti  ■  ij  : 


hence 


*0  -  - 


>*i2 .  r.  i-**.. 

— — -  +  U  ti  +  -  V 

>  1  K 


n  (1  -  Ktj) 


P  1  M1  1  2  U  f  M,  M,  M  " 

=  ~  T  k  (1  ~  ^  +  k  ^  (1  -  s;5  -  s;  n  mT  (3"] 

Taking  now  the  moment  of  burn-out  inturn  as  t  =  0,  the  velocity  at 
any  time  t  after  the  burn-out  is  given  by  the  expression 


V  =  V  -  g  t 
max  yo 


The  maximum  height  is  attained  at  the  time  t2  after  the  burn-out 


t  -  Vp|ax  -  1  n  Ml\  +  u  M° 

*2  •  —  -  -  i  (1~K ’  5„  ' "  B7  - 

The  value  of  this  height,  beyond  the  burn-out  point  is 

V2 

*  max 

A  *2  =  T“  * 


(3-14) 


Substitute  the  value  of  V„ 


B  [  u  .  !  (1  -  !i) 

2g0  [  M1  K  Mo 


(3-15) 


the  total  maximum  height  will  bes  7 0  +  &  Y2 
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Hence 


3.5  As  an  example  the  ascent  equations  have  been  applied  to 

the  computation  of  the  take-off  from  the  ground  level  for  a  lunar 
excursion  module  (LEM)  of  weight  4  tons.  Ail  weights  in  this 
problem  are  measured  at  the  standard  conditions  on  Earth. 

Assume  I  =  300  sec,  burning  time  =  40  sec,  a 

0 

constant  propellant  flow  W  =  c  =  20  lb/sec; 


u  =  I  ge  =  (300X32.174)  =  9652  ft/sec. 

c  _  1  1  W0  8000  WL  7200 

Hence  K  =  ^  =  7200  -  1.111,  ^  8000 


0.900. 


I 
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Bura-out  conditions 


y  =  V 
b.o  max 


.  ifi  ( 

K 


?1)  ♦  I 


t 


b.o. 


♦  I  ge  (2.306)  log1Q  Jjs 

=  -  (40X5.3)  +  ( 300)  (32.174H2. 3026 )  log^  (1.111) 
=  -  212  +  (9650X2. 3026X0. 04571)  =  -  212  +  1017 


=  805  ft/sec. 

Altitude  y.  attained  at  burn-out: 
b.o 


=  -  V-  + 1  s 

Jb .0  2  ue 


t  -  ( -  -  t )  In  — 
K 


(5.3) (1000) 


+  (9650) 


40  -  (400  -  40)  In  1.111 


=  -  4240  +  (9650) 


40  -  (360) ( . 1052) 


=  -  4240  +  (9650) (40  -  37.9)  =  -  4240  +  20260 


=  16020  ft. 

Maximum  acceleration  attained  at  burn-out 


—  +  K^I  = 
9e  9e  W1 


300  W 

=  -  0.164  +  -  0.164  +  (0.75X1.111) 

400  Wj 

=  0.670 

y  =  (0.670)02.174)  =  21.55  ft/sec2 
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The  time  of  coasting  to  the  maximum  altitude 


Vmax  806 

At  =  -2LS2  =  -  =152  sec. 

g0  5.3 

total  time  of  ascent  t  =  t  +  At  =  40  +  152  =  192  sec. 

max  b.o 


extra  height  beyond  the  b.o. 


Ay  = 


'max 


(005)' 


648025 


=  61100  ft. 


total  y  =  y. 

max  b.o 


2g  10.6  10.6 

+  Ay  =  16020  +  61100  =  77120  ft. 


This  result  is  checked  by  using  the  equation  (3-16)  which  now 
assumes  the  form 

,  w  2 

W  =  (8-78)  x  10  Un  wf*  + 


+  (3.86)  x  106  (1  -  -  -tn  ^  ). 

W0  Wi 

2  2 

y  =  ,-2.6.5?} .  (.1052)  +  (9650)  ( 400)  ( 1  -  0.9-0.1052) 

max  10.6 

(9.312)  x  107  x  1.1067 
10.6  x  102 

-  (9.65)  x  (4)  x  10^  x  (5.2)  x  10~8 
=  9.72  x  104  -  2.06  x  104  =  76600  ft. 

w 

The  same  pattern  is  applied  to  other  values  of  and 

"1 

the  results  are  presented  in  the  Table  IV.  All  computations  have 
been  performed  with  slide  rule  accuracy,  this  explains  the  small 
discrepancy  in  the  values  of  ymax  • 


I 
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CHAPTER  IV 


Vertical  Descent  and  Landing  Constant  Thrust 


4.1  The  problem  is  not  a  simple  reversal  of  the  take-off  problem, 
but  presents  some  different  aspects. 

Assume  that  at  a  certain  altitude  y0  the  horizontal  velocity 

of  the  vehicle  has  been  killed,  but  it  still  possesses  a  vertical 
velocity  y0.  This  velocity  plus  the  velocity  acquired  in  the  free 
fall  from  the  altitude  y0  must  be  reduced  to  zero  at  the  altitude 
just  above  the  ground. 

A  particular  case  is  when  y0  =  0.  It  is  shown  in  the  next 

Chapter  that  in  order  to  minimize  the  propellant  comsumption  it 
would  be  advisable  to  let  the  vehicle  to  drop  freely  for  considerable 
time,  and  then  to  apply  one  or  few  bursts  of  thrust.  This,  however, 
implies  a  re-ignition  of  the  engine  and  the  safety  considerations 
require  the  avoiding  of  such  procedure,  as  the  failure  of  the 
engine  to  ignite  would  be  fatal,  this  consideration  has  precedence 
over  the  fuel  economy.  Hence  we  examine  first  the  descent  under 
continuous  thrust,  which  may  or  may  not  be  modulated. 

Consider  first  the  case  of  constant  thrust  and  a  constant 

•  # 
propellant  flow  W  =  K  W0 .  Initial  conditions  are  y0  and  y0  =  V„. 

Final  conditions  ares  yf  =  0,  y  =  0.  More  properly  y^  is  a  small 

altitude  above  the  ground,  at  which  the  vehicle  may  hover,  giving 
to  the  pilot  the  opportunity  to  select  another  landing  place  if 
necessary.  At  any  time  t  since  the  beginning  of  the  maneuver  we 
have 


•  • 


y 


I  K 

go  + 


1  -  Kt 


(1) 
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First  and  second  integrations  yield 


y 


y 


vo  -  00  1  +  1  9e  — -  *  (2) 

1  -  Kt 

v»t  -  V- + 1  ’•[*  -  *> u  rnr]  + 


+  y0  <3) 

The  substitution  of  final  conditions  into  the  equations  (2)  and 
(3)  yields  two  equations  which  theoretically  can  be  solved  to  find 
K  and  t,  that  is  the  rate  of  fuel  consumption  and  the  time  of 
descent . 

Let  In  ■  *  .  =  X  , 

1  -  Kt 


V0  “  90  1  +  1  9e  x  =  0  > 

2 


(4) 


V0  t  - 


9o  1 


+  1  9, 


t  -  -  t)X 

K 


y0  =  o 


Note  that  X  >  0 


From  (4)  X  =  - 

Substitute  into  (5) 


^o  —  9 o  f 
1  9e 


(5) 


(6) 


v„  .  -  !g-l  +  y0  *  I  9e  t  ♦  (i  -  t)(T„ 


Solving  (6)  and  (7)  with  respect  to  K 

h2 

(I  9e  +  ^0)t  “  Fo 


90  t} 
(7) 


=  0 


So1 


i  2 

K  ” 


V  -  g  t 
o  ao 
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1 

K 


9e  1  +  y0 


3o 


t 


V 


o 


K  = 


g  t  -  v 
a  o  o 


=  f(t) 


(8) 


g  t* 

s  c 


+  i  gP  1  +  y. 


and  from  (6) 


g  t  -  v 

o  o 

_i _  e  i  ge 


1  -  Kt 


g„t  -  V 
■*  o  o 

Plotting  curves  K  =  f(t)  and  K  =  $(t)  their  intersections  yields 
the  values  for  t  and  K.  The  method  is  illustrated  by  considering 
*  vehicle  descending  from  the  altitudes  y0  =  5,000  ft,  10,000  ft, 
20,000  and  40,000  ft,  to  the  lunar  surface.  In  all  cases  it  has 
been  assumed  V0  =  -  400  ft/sec  and  I  =  300  sec.  The  curves  for 

the  solution  of  equation  (7)  have  been  computed  in  Tables  V, 

VI,  VII  and  VIII,  and  the  results  presented  on  the  Graph  A,  where 
the  solutions  have  been  obtained- 
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From  the  graph 
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l 


(10) 


s  ec 
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The  Tables  V,  VI,  VII  and  VIII  provide  data  for  graphical 
solutions  of  equations  (4)  and  (3).  For  all  altitudes  one  common 
curve  K  =  tj>(t)  is  plotted,  which  does  not  depend  on  y0.  Then  for 

any  value  of  y0  a  curve  K  =  f(t)  shown  by  equation  (8)  is  platted. 
The  intersection  of  i]i(t)  and  f(t)  yields  the  value  of  K  and  that 
of  tj  -  time  of  descent.  Some  obtained  results  have  been  checked 
by  substitution  into  the  equations. 


(1) 

t 

134 

mm 

98 

sec . 

(2) 

K 

8.2 

8.93 

9.4 

x  10-4  — i- 
sec 

(3) 

Kt 

0. 109 

0.0992 

0.0921 

(4) 

1  -  Kt 

0.891 

0.9008 

0.9079 

(5) 

1 

1  -  Kt 

1.1223 

1.1101 

1.1014 

(6) 

10910  1  -  Kt 

0.04999 

.04532 

.04198 

(7) 

( 2 . 3026 ) (6 ) 

0.1151 

0.1043 

0.0967 

(8) 

%  1 

710 

588 

520 

ft/ sec 

(9) 

9o  1  -  Vo 

1110 

988 

920 

ft/  sec 

(10) 

(9) 

i  ge 

0.1150 

0.1024 

0.954 

This  table  served  to  check  the  graphical  solutions  by  substitution 
into  the  equat i on  (6) wri tt en  as 


i  g  t  =  V 

— L_  =  12— - 

1  -  Kt  I  g 

.7 

Line  (^)  yields  the  values  for  the  left  side  and  the  line  (10)  for 
the  right  side  of  this  equation. 
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Cheek  of  the  equation  (5),  written  as 


1  -  Kt  ,  x 
t  -  — — —  in  ( ) 
K  l  “  Kt 

we  obtain  for  the  left  side  in  the  case  y0  =  lOOOO  ft. 

(5.3)  (12321) 

. .  -  .  -  10,000  +  400  x  111  = 

2 

=  32630  -  10,000  +  44400  =  67030  ft; 


So 


“  To  “  V 


=  I  Se 


and  for  the  right  side: 


I  9< 


=  9652 


1  -  Kt 

t  -  In  ( 


1 

1  -  Kt* 


0.891  x  104  ( 

111 - -  ... - —  (0.1043) 

8.93 


=  9652 


111 


=  9652  (111 


8-91  3  , 

-  x  10J  (0.1043) 

8.93 

104.07)  = 


=  (9652)  (6.93)  =  68200  ft. 


The  agreement  is  fairly  good  for  the  slide  rule  accuracy.  The 

results  of  the  computations  are  summarized  in  the  following  table 

and  graph,  where  t  is  the  time  of  descent  from  the  initial 

W0  “  W, 

altitude  y0  to  y  =  0,  Kt  =  — n -  is  the  ratio  of  the  propellant 

"o 

consumed  during  the  maneuver  to  the  initial  weight  of  the  vehicle 
at  y  =  y0.  The  magnitude  of  the  initial  and  final  accelerations 
are  computed  from  the  equation 


0.164  +  KI 
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•  * 

9e 


W, 


» 


where 

V1  =  W0  (1  -  Kt)  f 

for  t  =  0,  —  =  -  0.164  +  KI  =  initial  acceleration  in  terms  of  the 
9e 


Earth  standard  gravity  ge.  , 

if  ° 

For  each  landing  case  theVfuel  consumption  and  the  thrust 

are  constant.  This  implies  that  whatever  is  the  fuel  flow,  the 

specific  thrust  I  remains  constant.  The  ratio  of  thrust  to  the 

F  W 

initial  weight  i s  ~  I  =  KI .  The  weight  here  is  the  standard 

»o  wo 

Earth  weight  of  the  vehicle. 

The  results  obtained  for  descent  from  different  altitudes  are 
summarized  in  the  Table  IX  and  presented  graphically  on  the 
Graph  B. 
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Note  that  the  thrust  F0  is  compared  to  the  initial  weight  of 

the  vehicle  at  the  Earth's  standard  gravity.  If  the  lunar  initial 
weight  is  considered,  then 


(wA) 
o  e 


9e  Fo  1 
—  =  (-2)  - 


("oh  g0 


W0  0.164 


We  have  investigated  the  case  of  a  descent  from  different  altitudes, 
with  the  same  initial  vertical  velocity.  Consider  now  the  case  of 
a  vertical  descent  from  the  same  altitude  but  with  different 
initial  velocities  VQ.  For  an  example  the  altitude  y0  =  10,000  ft 
has  been  selected,  and  for  V0  the  downward  values  0,  -400,  -800, 
-1000  ft/sec  assumed.  The  equations  (8)  and  (9)  have  been  used  for 
plotting  the  curves  K  =  f(t)  and  K  =  \j»(t);  the  coordinates  for 
these  curves  are  computed  in  the  Tables  X,  XI,  XII  and  V,  the 
curves  are  plotted  on  Graphs  C  and  D,  yielding  graphical  solutions. 
The  results  are  summarized  on  the  Table  XIII;  then  K,  Kt,  and 

% 

on  the  Graph  E,  which  shows  that  for  each  initial  altitude  there  is 
an  initial  V0  yielding  a  minimum  propellant  consumption 


final  acceleration  are  plotted  versus  the  initial  velocity  V( 


v.  propellant  weight  *  .  .  .  ,  .  . 

Kt  =  *■—£«: — — - —  ” .  In  the  considered  example  this  minimum 

initial  weight 

occurs  for  y0  =  10,000  ft,  V0  =  -  570  ft/sec. 
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I 


360  sec.  K  =  4.99  x  10~4  =1=  Kt  =  0.1798 

sec 
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CHAPTER  V 

Intermittent  and  Modulated  Thrust 

5.1  Assume  that  the  horizontal  velocity  of  the  vehicle  has  been 

killed  at  the  altitude  y0.  Its  vertical  velocity  is  y  =  V0,  which 

may  be  zero  in  particular  cases.  In  the  previous  chapter  a 
vertical  descent  under  a  constant  continuous  thrust  has  been 
considered,  without  the  stalling  and  restarting  the  engine.  Assume 
now  that  the  engine  has  been  switched  off,  the  vehicle  in  free  fall 
will  reach  the  altitude 


yi  =  y0  +  vo  t 


g0  t2 
2 


> 


with  a  vertical  velocity  yj  =  V0  -  gQ  t.  At  this  altitude  the 

engine  will  be  again  ignited,  then  the  rate  of  the  propellant  flow 
will  be  computed  as  it  has  been  shown  in  the  Chapter  IV  in  order 
to  reach  the  ground  with  zero  velocity. 

Without  embarking  on  a  generalized  optimization  problem  let 
us  compare  two  methods  of  descent  from  the  point  of  view  of 
propellant  economy.  Consider  a  descent  from  the  altitudes 
y0  =  20,000  ft  and  y0  =  40,000  ft,  assuming  the  initial  vertical 


V0  =  -  400  ft/sec  for  both  cases.  Assume  that  the  ignition  occurs 
when  the  vehicle  reaches  the  altitude  yi  =  10,000  ft.  Graph  on 


Figure  E  enables  us  to  find  K  and  fuel  consumption 


Wq  -  Wj 
Wo 


=  Kt  required  to  reach  the  ground. 
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Powered  descent  from  20,000  ft:  (as  given  in  the  Table  IX) 

Propellant  consumption  Kt  =  0.109  of  the  initial  weight 
Time  of  descent  t  =  134  sec. 

Maximum  acceleration  =  .112  ge 

Drop  from  y0  =  20,000  ft  to  =  10,000  ft.,  time  of  the  fall  is 
computed  from  the  equation 

g0  t2 

Ay  = - - -  +  V0  t, 

where 

Ay  =  yx  -  y0  =  -  10,000  ft. 


Hence 


gt  -  2  V„  t  +  2  Ay  =  0; 


t  =  At  = 


-  2g  Ay 


Considering  only  the  value  At  >  0,  we  obtain 

+  y40^  +  2x  5T3X  10^000 


At  = 


-  400 


5.3 


116.3 

5.3 


The  velocity  attained  is 


Vl 


2g0  Ay  =  -  516.3  ft/sec. 


21.95  sec. 
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In  the  subsequent  powered  descent  we  obtain  from  the  Graph  E 
based  on  Table  XIII. 

t  =  81  sec.  ,  t  +  At  =  103  sec. 


Kt  =  0.0935 


K  =  12.1  x  10 


-4 


sec 


Max.  accel.  =  0.24  ge 


Both  methods  are  compared  in  the  following  table: 


Table  XIV 


Descent  from  y0  =  20,000  ft.,  VQ  =  -  400  ft/sec. 
_ I  =  300  sec. _ 


Constant 

Thrust 

Intermittent 

Thrust 

t  Total 
time  of 
descent 

134 

103 

sec 

"o 

0.109 

0.0935 

y/g 

.112 

.24 

K 

CM 

• 

CO 

12.1 

x  10"4  — 
sec 

Hence  the  second  method  shows  a  not  negligible  fuel  economy,  with- 
an  acceleration  well  below  the  Earth's  gravity. 


109  -  93.5  =  16.5 
109 


15.1  p.c. 


fue  1 
economy 


109 
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The  similar  computation  applied  to  a  powered  descent  from  the 
altitude  40,000  ft.  yields: 

At  =  56.8  sec.  Vj  =  691.4  ft/sec. 


Table  XV 


Descent  from  y0  =  40,000  ft.,  V0  =  -  400  ft/sec. 
1  =  300  sec. 


Constant 

Thrust 

Intermittent 

Thrust 

t 

168 

56.8  +  50  =  106.8 

sec . 

*1  > 

o  1* 

II 

.125 

.096 

Y/Qe 

.092 

0.54 

K 

7.44 

20.8 

x  10-4  -i- 
sec 

In  this  case  the  economy  of  propellant  is  even  better,  being 


.125  -  .096 
.125 


•  £  7 

-  =  23.2  p.c. 

125 


The  max.  acceleration  is  still  only  about  half  of  the  Earth's 
gravity. 


5.2  Modulated  Thrust 

In  order  to  avoid  the  stalling  of  the  engine  and  to  obtain 
some  fuel  economy  in  conditions  similar  to  those  with  the 
intermittent  thrust,  a  modulated  thrust  can  be  used.  Modulation 


I 
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is  obtained  by  varying  the  flow  of  propellant.  This  involves 
extra  difficulties,  as  this  flow  should  be  continuously 
nodulated.  Besides  we  assume  a  constant  specific  thrust  which 
is  the  measure  of  conversion  of  the  thermal  energy  released  in 
a  chemical  reaction  into  the  Kinetic  energy  of  the  mass  flow. 

This  conversion  is  optimum  for  a  certain  value  of  the  mass 
flow  corresponding  to  a  given  size  and  shape  of  the  nozzle  and 
combustion  chamber.  The  value  of  I  may  drop  if  these  conditions 
are  not  fulfilled,  hence  an  adjustable  nozzle  will  be 
required.  This  entails  extra  complexity  of  design,  which  should 
be  rather  avoided. 

So  we  have  either  to  assume  that  I  =  const.,  or  provide 
the  diagram  of  its  variation  with  the  mass  flow.  In  this 

section  we  assume  I  =  const  =  300  sec.,  irrespectively  of  the 

» 

magnitude  of  W  or  of  the  size  of  the  nozzle. 

The  variation  of  the  thrust  F  is  assumed  to  follow  the  law 
F  F  Q 

-  =  —  =  const;  where  W  is  the  weight  of  the  vehicle  at  a  given 

W  Wq 

time  t.  Since  the  gravity  acceleration  is  assumed  to  be 
constant,  the  decrease  of  the  weight  is  due  only  to  the  propellant 
flow.  W0  is  the  weight  of  the  vehicle  at  the  beginning  of  the 
maneuver,  that  i s  at  y  =  y0.  Its  value  is  the  initial  weight 
of  the  vehicle  measured  on  the  lunar  surface,  (not  on  the  Earth). 
For  a  given  altitude  y0  and  vertical  velocity  V0  we  have  to  find 

9 

w 

the  time  of  descent  tj  ,  the  modulation  of  K  =  —  as  a  function 

"0 

of  time  and  the  ratio  of  the  propellant  consumed  to  the  initial 
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weight  of  the  vehicle 


AW  fll 

!  rr  I. 


Kdt.  We  start  with  the  equation 


of  notion 


W  •* 

—  y  =  F  -  W 
9o 


or 


*  *  I  «  ‘  1 


_  JFo 

9o  “  W^  "  I  9o  =  coast* 


The  first  and  second  integrations  yield 


(5-1) 


y  =  0o 


Fo  |  •  F0 

_  .  1  t  +  yo  =  -  -  !  ,  go  t  +  V0 


y  =[  L  -  l\  +  V  t  +  y 

y  I  W  2  0  y< 


(5-2) 


(5-3) 


At  the  ground  level  we  ought  to  have 

y  =  0,  y  =  0  and  obtain  two  equations  with  the 

unknowns 

J  F0 

ti  and  -»  . 


From  (5-2) 


ti  =  - 


r-  -  1  9o 


The  (5-3)  yields 


hence 


and 


21  -  -  1 

W„ 


—  -  1 

W„ 


F  V  2 
r  o  _  v  o 


+  ya  =  o 


+  i 


-  i  = 


ti  = 


wo  2  y0  g0 

2yo 

V0 


V  2 
_  *0 _ 


2  Vo  00 
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The  rate  K  =  —  of  the  fuel  flow  can  be  found  from  the  equation 
"o 


W 

W  =  Wo  -  J  -B-  dt  ,  K 

"o  "0 


is  variable  in  this  case 


.  ge  .  dW 

F  =  IW  —  ;  but  W  =  -  — 

90  dt 


g~  dW 

F  =  -  I  —  — -  ;  divide  both  sides  by  W. 

g0  dt 


F  _  [o  _  g*  «  1 
w  =  w0  =  “  g0  '  dt  *  w 


dW  =  _  go  1  Fo  dt 

W  ge  i  w0 

for  any  given  time  t  which  elapsed  since  the  beginning  of  the 
maneuver  we  have 

•  Wo  go  1  Fo 
"  W  =  ge  i  W0 

.  9o  1  Fo 

ana  —  ▼  fiP  1 

W  =  W0  e“  9e  0 

The  amount  of  the  propellant  consumed  since  t  =  0  till  t  =  t 
is  90  Fo  x 

«F  _  «0  -  W  ~  ~e  v~0  i 

—  —  ■  -  i  -  e 


®o  Fo  1 

K  =  -l  "  =  + 

Wo  dt  ge  W0  I 


l 
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As  an  example  we  compute  the  descent  from  y0  =  40,000  ft, 
V0  =  -  400  ft/sec^for  I  =  300  sec  and  compare  the  results  with 
those  for  constant  and  intermittent  thrust. 

2  y  2  x  4  x  104 

t .  =  =  — —  =  200  sec. 

V0  4  x  102 


F0  V02  4  x  4  x  104 

—  =  — -  +1=1+  - 

wo  2  y0  8o  2  x  4  x  104  x  5.3 


4 

=  -  +  1  =  1.377 

10.6 


The  final  value  of  the  variable  thrust 

£o  Fo  li 

.  Jj  !i » £*  !i  1  r 

w0  Wj  w0  w0  w0  w0 


—  —  -  «  (0. 164)0. 377)(  )  =  7.53  x  lo“4  -i.. 

ge  W0  I  300  8ec 

for  tj  =  200  sec,  we  have 


Hence 


-  7.53  x  10-4  x  2  x  102  -  1.506  x  10"1 

e  =  e 


-0. 1506 

=  e 


0.8602 


F1 

—  =  ( 1.377) ( .8602 )  ■  1.184 
"o 


/ 
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Propellant  consumed  in  the  descent 


e-0.1506  _  j  _  0#8602 


.1398; 


-0.1506 

e 


1.1625 


Acceleration  is  constant  in  this  case: 


0.164 


+  !l2  22 

So  5e 


I  = 


=  -  0.164  +  0.226  =  0.062. 


TABLE  XVI 

K  =  7.53  x  10“4  x  e“7*53  *  10_4  1 


t 

0 

50 

100 

150 

200 

sec 

7.53  t  x  10~4 

0 

0.03765 

0.0753 

0.1130 

0.1506 

e-7 .53  x  10“4t 

1 

0.9630 

.9275 

.8932 

.8602 

K=  L 

W0 

7.53 

7.25 

6.99 

6.73 

6.48 

—4  1 

X  10  9  -i- 
sec 

The  table  provides  the  programming  of  the  fuel  flow  during 
the  descent  from  y  =  40,000  ft.  It  is  valid  only  if  the  specific 
thrust  I  =  300  sec.  remains  constant  in  spite  of  the  varying 


conditions  in  the  nozzle  and  combustion  chamber.  If,  however,  I 
is  not  constant,  its  variable  values  provided  by  the  engine  data 
must  be  introduced  in  the  expressions  for  thrust  and  rate  of 
propellant  flow. 

The  results  of  the  present  investigation  can  be  summarized 


in  the  following  table 
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Hence  a  modulated  thrust  did  not  produce  an  expected  reduction 
in  propellant  consumption,  because  the  lower  rate  of  fuel  flow 
was  offset  by  the  length  of  the  descent  time.  Besldes^this  case 
t  would  require  a  continually  adjustable  nozzle  in  order  to  maintain 

a  constant  specific  thrust.  This  would  entail  extra  complexity 
of  the  design  and  a  decrease  in  reliability.  So,  at  least  for 
the  given  initial  conditions,  the  choice  would  be  rather  between 
the  cases  A  and  B. 


